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Abstract. The explicit energy dependence of the single particle self-energy (dispersive effects), 
due to short range correlations, is included in the treatment of neutron matter superfluidity. 
The method can be applied in general to strong interacting fermion systems, and it is expected 
to be valid whenever the pairing gap is substantially smaller than the Fermi kinetic energy. The 
results for neutron matter show that dispersive effects are strong in the density region near the 
gap closure. 
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1. Introduction 

Neutron and nuclear matter superfluidity is one of the main issue in the physics of 
neutron stars. Superfluidity is expected to play a major role in some of the most strik- 
ing phenomena occurring in neutron stars, like glitches and post-glitches transients 
vortex pinning ||, neutron star cooling, and maybe strong magnetic field penetration 
H]. However, since the observational data are only indirectly related to superfluidity, or 
need explicit models for their interpretation, a firm theoretical prediction of the super- 
fluidity strength, based on microscopic ab initio calculations, appears highly required. 
Unfortunately, neutron and nuclear matter are strongly correlated systems, where short 
range correlations dominate the overall interaction energy, even at densities well below 
the saturation value. The superfluidity problem turns out, therefore, to be a complex 
many-body problem, where the delicate balance between short range interactions and the 
long range pairing correlations needs an accurate treatment. This problem was firstly 
considered in the works of ref. Q, within the variational Jastrow method. The medium 
effect on the effective pairing interaction was investigated in refs. [||], |6|. The correlation 
effects were treated in ref. [0] within a generalization of the Babu-Brown approach || to 
the effective nucleon-nucleon (NN) interaction and in the weak coupling limit. In general, 
these microscopic approaches seem to indicate a reduction of the pairing gap due to the 
medium, with respect to BCS approximation with the bare interaction. The use U in 
the BCS scheme of realistic bare nucleon-nucleon interactions, which reproduce the ex- 
perimental phase shifts, can be a good starting point for a more sophisticated many-body 
treatments, and the connection between the pairing gap value and the phase shifts has 
been elucidated, in general, in ref. ]T0| . 



In all these microscopic approaches the single particle spectrum is usually considered 
within the effective mass approximation, or taken from normal state calculations. Dis- 
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persive effects, due to the energy dependence of the single particle self-energy, are usually 
neglected, or considered only in the weak coupling limit, and no first principle scheme to 



include them in the general gap equation has been proposed. Only recently [1 1| a self- 
consistent scheme has been developed, where the short range correlations and the pairing 
problem are treated on the same footing. The method is numerically complex, and it has 
been solved only for a schematic interaction ||11|| . It would be desirable to understand in 
simple terms the effect of the self-energy dispersion on the pairing strength, also in the 
case of strong coupling as in neutron matter. In this letter we present a general scheme 
for including dispersive effects in the gap equation, which is simple and accurate, provided 
the gap is substantially smaller than the Fermi kinetic energy. We then apply the method 
to neutron matter superfluidity with realistic interaction, and we show numerically that 
the size of the effect can be large in the vicinity of the gap closure. 

2. Including self-energy in the gap equation 

The general many-body theory of pairing in fermion systems has been formulated just after 



the BCS [0 solution has been introduced, and it can be found in standard textbooks ||T3 
14]. The main framework is the Green' s function (GF) formalism, which generalizes the 
Gorkov' s method beyond the BCS approximation. The single particle Green's function 
Q has a 2x2 matrix structure, with normal diagonal components Fi and abnormal off- 
diagonal components F 2 , 

G(k ,A- ( F ^ k >u) Hb") \ c-Uh,A-( h~u + M(k,uj) A(k,u) \ 
y{K ' U) \ F 2 (k,uj) - Fi(k, -oj) J ' [K ' UJ) { A{k,u) -(i h + u + M{k,-u)) J 

(1) 

In the expression for the inverse Green's function Q l , we have introduced the quantity 
€k = h 2 k 2 /2m — n as the single particle kinetic energy, with respect to the chemical 
potential fx, the diagonal single particle self-energy M(k, u) and the momentum and 
energy dependent gap function A(k, to). Here we are assuming s-wave singlet pairing, and 
therefore we omit spin indices. They simply express the coupling between the time-reversal 
states {k, |) and (—k, j). Both M(k, to) and A{k, u) can be expanded in terms of the NN 
interaction and the full GF itself, which in general entails a self-consistent procedure. The 
gap function A{k,uj), however, is solution of the homogeneous generalized Bethe-Salpeter 
equation [TBI], and as such it satisfies the generalized gap equation 



r , I(kuo,k'u')A(k',u') 

where I (fad, k'u') is the irreducible NN interaction at zero total energy and momentum. If 
one takes the bare NN interaction for the interaction I, and the Hartree-Fock approxima- 
tion for the diagonal self-energy M(k,u), the standard BCS approximation is recovered. 
It has to be noticed that the energy dependence of the gap function A(k,u) originates 
only from the energy dependence of the irreducible interaction I. In fact, if the interac- 
tion is taken as energy- independent, the gap function is also energy- independent, despite 
the possible energy dependence of the self-energy M(k,u). Since we are looking for an 
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estimate of the dispersive effects, we indeed assume the irreducible interaction as energy 
independent, while we keep the full energy dependence of the self-energies. Then, the 
energy integration appearing in the gap equation (Q) can be performed with a good ac- 
curacy in the limit of small self-energy imaginary part, since than the main contribution 
is expected to come from the poles close to the real axis. The denominator is an even 
function of the energy u (we remind again that single particle energies are measured with 
respect to //), and, therefore, the kernel presents two poles, symmetrical with respect to 
the origin in the complex tu-plane. This is a feature typical of the superconducting phase. 
Formally, the pole energies ±E k are the solutions of the implicit equation 



±E k = -{M(k, ±E k ) - M{k, TE k )) ± ■ 



i 2 



e k + -(M(k,-E k ) + M(k,E k )) 



+ A(ky (3) 



If the energy dependence of M(k, uj) is neglected, than Eq. (|3]) reduces to the usual square 
root expression for the quasi-particle excitation energy of the BCS approximation. On 
the other hand, in the non-superconducting limit A — > 0, and neglecting the imaginary 
part of M(k, uj), one can verify that Eq. (||) reduces to the usual self-consistent equation, 
e.g. Brueckner |T5|], for the single particle energy e k 



e k = e k + M(k,e k ) 



(4) 



Equation is valid whenever A is negligible, in particular for momenta far away from 
the Fermi surface, since then \A k \ <C |efc|. We now make use of the assumption of a 
gap A smaller than the normal self-energy, which is mainly determined by short range 
correlations. The size of the normal self-energy is indeed of the same order of the Fermi 
kinetic energy Ep. In this case, Eq. ([|) will be valid to order A{kp) / Ep, and therefore 
on the right hand side of Eq. (H) we can replace E k with e k , solution of Eq. (ffl), to get 



Eh 



-(M(k,\e k \)-M(k,-\e k \)) + 
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(M(k, —e k ) + M(k, e k )) 



A(k) 2 (5) 



where the self-energy is now calculated in the normal phase. The procedure is justified, 
provided M(k, uj) is a smooth function of uj. It will be further discussed below. 
The residue R of the kernel at each one of the pole can be easily calculated 



R 



a k 




i--((i-e k )a k + (i + e k )b k 

_ e k + k(M(k,-E k )+M(k,E k )) 



E, 



w=E k 



'dM s 
duj 



(6) 
(7) 
(8) 



In the limit A — > 0, Eq. (6) is the usual expression for the quasi-particle strength 
Z k) provided the momentum k is close enough to kp. The corrections to the normal 
phase value of 1 — R^ 1 are of the order A(kp)/Ep, and therefore the residue R can be 
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identified with Z^, at least in the vicinity of the Fermi momentum. Far away from the 
Fermi momentum, the residue R has still the expression of Eq. (6) in the limit of small 
imaginary part, despite the quasi-particle concept becomes less meaningful, since its width 
can be large (but it can be still much smaller than the real part of the energy). In this 
case the procedure is just an accurate method of calculating the energy integral (pole 
approximation). For simplicity, the residue will be denoted by Z k in all cases. Within 
these approximations, the generalized gap equation @ reads 

A(*) = - £ I(k, k')Z k , Ak ' (9) 



sj [h + \(M(k, -e k ) + M(k } e k ))] 2 + A(k) 2 



k ' 2 1 

Eqs. (|), (6), (H) contain the main result of the paper. It has to noticed that in the 
generalized gap equation (Q) the square root in the denominator does not coincide with 
the quasi-particle energy of Eq. ([5]), in contrast with the usual BCS approximation. This 
feature is general and it is not bound to the approximation of Eq. (|5[). The approximation 
of Eq. (|) to the expression inside the square root is valid up to corrections of order A/E F . 
They can be absorbed in large part by a small shift of the chemical potential fx. 

Before going to the application of the formalism to neutron matter, let us discuss Eq. 
(§) in the extreme weak coupling limit, where one assumes that the main contribution to 
the momentum integral is concentrated around the Fermi surface. In this limit, following 
the standard procedure of expanding the integrand of Eq. (||) around k F , one gets |16| 



F m F ^ m F Z F 7r 2 n I(k F )^ ^ 

where no is the density of state for the free Fermi gas and m F the so called k — mass 
(in units of the bare mass) [[L7j]. The interaction I(k F ) is the diagonal matrix element 
of the NN potential in the considered channel ( 1 S'o for neutron matter), in the plane 
wave representation. The self-energy effects are, therefore, contained mainly in the factor 
m F Z F , which can be written also as m*Z F , since the full effective mass m* = m F /Z F 



17|j . This is the standard result for the weak coupling limit [[Tj|. Eq. (H) generalizes 
the treatment to the case where the contribution from momenta far from the Fermi 
momentum is relevant. The appearance of the k — mass is a peculiar feature of the 
pairing phenomenon, and it is a direct consequence of the coupling between time-reversal 
states. In Eq. (|9[) the combination M(k, —u) + M(k, u) gives rise to the combined density 
of state of the pair {(k,u); (—k, —u})}, which is mainly determined by the k — mass. 



The weak coupling limit is not valid in general for neutron or nuclear matter [T^j, if 
one starts from the bare NN interaction. This can be seen directly from the observation 
that often the gap equation has a well defined solution even when the interaction matrix 
element I(k F ) is positive. This is due to the dominant role of the off-diagonal matrix 
elements I(k,k'). Therefore, one must solve the more general equation (0) in this case. 
The above considerations are, anyhow, still valid. 
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3. Application to neutron matter superfluidity 



In order to estimate the dispersive effects on the superfiuid gap of neutron matter, we 
have solved Eq. (^|), with the bare Argonne V14 potential as the pairing interaction I(k, k') 
and with the self-energy calculated in the Brueckner approximation at the lowest order, 
see Fig. la, with the same interaction. The higher order contribution of Fig. lb turns out 
to be indeed negligible in the relevant density range. In the superfiuid phase, in principle, 
the diagonal self-energy M(k,u) differs from the self-energy in the normal phase. The 
main contribution not present in the normal phase originates from the coupling of the 
single particle motion with the superfiuid collective modes. The latter correspond mainly 
to the center of mass motion of the Cooper pairs and their possible "vibrations" Hl9| , pOfl . 
The branch starting at zero energy, in the long wave-length limit, is the branch of the 



Goldstone boson [19], corresponding to the gauge invariance symmetry breaking at the 



superfiuid phase transition. This contribution to the diagonal single particle self-energy 
is expected to be at most of the order of the superfiuid condensation energy per particle, 
and therefore negligible with respect to the typical short range correlation energy, as 
calculated e.g. in Brueckner theory, at least to the extent that A/Ep <C 1. For the same 
reason, the deviation of the occupation number from the free gas value and the presence 
of a forbidden energy region, of order A, around the Fermi energy, typical of the pairing 
phenomenon, seem to play no relevant role in determining the size of the self-energy. It 
appears, therefore, justified to adopt for M(k,uS) its normal phase value. 

The choice of the bare interaction for I(k, k') is suggested by the observation that no 
ladder summation should be included in the irreducible interaction kernel I(k, k') [TS|, |T5| ]. 
Of course, other terms, like polarization diagrams, should be included but here we 
want simply to single out the dispersive effects, and therefore it appears meaningful to 
compare the results obtained with and without self-energy, within the same scheme of 
approximation. 

In Fig. 2 is reported, for kp = 0.9/m -1 , the imaginary part of the neutron self-energy 
Mj(k,e(k)) at the quasi-particle pole, as a function of the momentum k, together with 
the real part of the quasi-particle energy e(k) (calculated with respect to /i). One can 
see that indeed the imaginary part is small with respect to the real part. The situation 
is completely similar for the other densities. The residue of Eq. (6), which appears 
in Eq. (^j) for the gap function, is reported in Fig. 3 for three densities. According 
to Migdal-Luttinger theorem ||18||, the value Zp of at k = kp is the discontinuity of 



the momentum distribution at the Fermi momentum (in the normal phase). One must 
have, therefore, < Zp < 1. One can see, however, that Z^ exceeds 1 slightly in some 
interval well above kp. This is not surprising, since for large u values, at fixed k, the real 
part of the self-energy is an increasing function of u, and it approaches asymptotically 
an energy- independent value (and therefore Zj. — > 1). As already mentioned, for 
large momentum the pole approximation is just an accurate method of calculating the 
relevant energy integral. The position of the pole is, of course, not exactly on the real 
axis, but numerical estimate of the second derivative of the self-energy shows that to 
calculate Z^ on the real axis is an extremely good approximation. The factor Zk is 
also related to the single particle occupation number n(k). In fact, the contour integral, 
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closed in the lower complex plane, of the single particle GF equals 1 — n(k), k > hp. 
Besides the pole, this integral receives contribution from the regular (non polar) part of 



the single particle GF [IS]. If Z(k) > 1, this means simply that the regular part leads to a 
slightly negative contribution. We have checked, indeed, that the single particle spectral 
function, calculated with the same self-energy (including of course the imaginary part), 
satisfies the sum rules and gives well defined occupation numbers n(k). A full account of 
the calculations will be reported elsewhere. Anyhow, this small deviation from 1 does not 
affect at all the final results, and one can take = 1 in this momentum region. 

Finally, in Fig. 4 is reported the pairing gap at the Fermi energy as a function of 
density for three different cases : i) without self-energy ||, ii) with the Z/~ factor in the 
numerator of the gap equation (|9]), iii) with both the self-energy in the denominator and 
the Zfc factor. The reduction of the pairing gap is substantial at the highest densities, 
near the gap closure. 

3. Discussion and conclusion 

We have developed a method to include the single particle self-energy in the gap equation, 
and in particular dispersive effects, beyond the usual BCS approximation. The method 
rely on the assumption of a small gap with respect to the Fermi kinetic energy and strong 
short range correlations, typical of the neutron matter in the inner crust of neutron stars. 
The results indicate that dispersive effects can strongly reduce the gap value near its 
closure. The effect is due both to the quasi-particle strength Z k and to the k — mass, 
which enter in the generalized gap equation (m. This result appears in line with the work 



of ref. [11], where the self-consistent treatment of pairing and short range correlations 



seems indeed to reduce strongly the gap value mainly because of these two factors [|Tl], ^1 
Of course, before drawing any conclusion on the pairing strength in neutron matter, 
one should include, along a consistent scheme, the correlation effects on the irreducible 
interaction I(k, k'). Work in this direction is in progress. 

The extension of the method to symmetric nuclear matter, possibly relevant for pairing 
in nuclei, appears problematic. In that case the approximation of a small imaginary part 



looks less justified away from the Fermi momentum ||22|| . Furthermore, the neutron- 
proton pairing in the 3 Si — 3 D\ channel is too strong, i.e. A ~ Ep, to be treatable in the 
proposed approximate method. In both self-consistent procedure, involving both 

the self-energy and the effective interaction, seems to be the only viable method. 



6 



Figure captions 



Fig. 1.- One (a) and two (b) hole-line diagrams contributing to the nucleon self-energy 
in the Bethe-Brueckner-Goldstone expansion. The wavy lines indicate Brueckner G- 
matrices. 

Fig. 2.- Real part e(k) and imaginary part Mj of the quasi-particle pole, as a function of 
the momentum k at the Fermi momentum kp = 0.9/m _1 . 

Fig. 3.- The residue Z k at the quasi-particle pole as a function of the momentum k at 
three values of the Fermi momentum kp. 

Fig. 4.- The superfluid gap value, at the Fermi momentum, as a function of density, in 
the case of free single particle spectrum (diamonds), with the inclusion of the factor 
(crosses) and with the inclusion of both Zk factor and the self-energy in the single particle 
spectrum (squares). 
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